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Abstract 

We give a simple criterion for checking the so called Edwards' hypothesis in 
certain zero-temperature, ferromagnetic spin-flip dynamics and use it to invalidate 
the hypothesis in various examples in dimension one and higher. 
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1 Introduction 

In many physical systems, the dynamics at low temperature or high density is so slow that 
the system is out of equilibrium at all practical time scales, so that ordinary thermody- 
namics does not apply to those situations. As a consequence, because of their practical as 
well as theoretical interest, non-standard thermodynamics have been proposed to describe 
out-of-equilibrium systems. (The literature concerning such systems is huge - for a few 
examples across the spectrum, see for instance [Ill21IIIlSII3IIE0IElin], an d references 
therein.) 

In the context of the slow compaction dynamics of granular materials, Edwards f^lTTl 
IT!?] proposed to compute thermodynamic quantities by means of a flat ensemble average 
over all the blocked configurations of grains with prescribed density, leading to a natural 
definition of configurational temperature. Since the approach is not justified from first 
principles, its validity has to be tested with specific models and experiments (see, for 
example, j7j and references therein). 

There are other situations in which the idea of an effective temperature is very ap- 
pealing, and it is tempting to extend the range of applicability of Edwards' idea to other 
systems with a large number of "metastable" states. In this broader context, what came 
to be known as Edwards' hypothesis consists in assuming that all the metastable states 
in which a system can be trapped are equivalent for the dynamics. This corresponds, as 
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before, to computing thermodynamic quantities using a flat ensemble average over all the 
metastable states. 

A situation where this prescription makes sense is represented by glassy dynamics, 
which are often described as a slow motion in a complex energy (or free energy) land- 
scape, with many "valleys" separated by barriers. Several approaches have been proposed 
to make this heuristic picture more precise, and different notions of "valleys" have ap- 
peared, but no general and unambiguous definition of metastable states has appeared yet. 
Nonetheless, once these valleys are appropriately defined, one can estimate their number 
Ne(N) at fixed energy (or free energy) density E for a fixed size N of the system. Ne^N) 
generally grows exponentially with the size of the system: 

tf E (N)~exp{NS B (E)), (1) 

where S C (E) is a configurational entropy or complexity. The key question concerning the 
dynamics is whether all valleys play the same role, or whether they have different dynami- 
cal weights, according to their basin of attraction. This question arises, for instance, when 
a system is instantaneously quenched into the glassy (low temperature) phase, starting 
from a disordered (high temperature) configuration. 

In this context, Edwards' hypothesis is valid for some mean-field models, where valleys 
are know to be explored with a flat measure ^3]. One can therefore define a thermody- 
namics based on the flat ensemble over valleys and a configurational temperature T c by 

1 = — • (2) 
T c dE y } 

Besides the mean-field case, another situation where valleys can be unambiguously 
defined is the zero-temperature limit, where no energy barrier can be crossed and the 
valleys correspond to the blocked configurations under the chosen dynamics. Models 
with constrained dynamics have been extensively studied in one dimension (see jTH] and 
references therein), where Edwards' hypothesis has been tested using exact results on the 
statistics of blocked configurations, as well as numerical simulations. The results obtained 
so far seem to imply that the hypothesis cannot be applied to zero-temperature spin-flip 
dynamics |15j . 

In this paper, we introduce a general criterion for checking Edwards' hypothesis for 
attractive (this refers to the ferromagnetic nature of the interaction between spins - see 
Section fTTTj) symmetric Ising dynamics with initial configuration chosen from a symmetric 
distribution that satisfies the FKG inequality (see II. II and [ini[T7J[TH]). As a first applica- 
tion (Section 12. lj) . we rigorously invalidate Edwards' hypothesis for the one-dimensional 
constrained Glauber dynamics analyzed in [TH] (see also [11112011211), which we use as a 
prototype. Later (Section l3.1|) . we show how the criterion can be easily applied to other 
dynamics in higher dimensions, where the analytic techniques of ^3] cannot be used and 
exact results are not available. 
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1.1 FKG Inequality and Harris' Theorem 



We give here the tools needed to check Edwards' hypothesis for attractive (see the next 
paragraph), zero-temperature Ising dynamics. In the context of an Ising spin model on 
a lattice L, we will call increasing an event £ such that its indicator function Is(cr) is 
increasing in the number of plus spins present in the configuration a = {cr x }x€iL, o x = ±1. 
If £\ and £ 2 are two increasing events, the FKG inequality (see, for instance, [T7[ ITHj ) 
states that, roughly speaking, the occurrence of £2 makes £\ more likely, or more precisely, 
the conditional probability of £\ given £ 2 is larger than or equal to the probability of £\. 



Many interesting distributions satisfy the FKG inequality (jHJ), among them are Gibbs 
measures (with some restrictions) and product measures, and in particular the symmetric 
Bernoulli product measure from which the initial configuration of the constrained Glauber 
dynamics of Section |2] is chosen (corresponding to a spin system prepared at "infinite" 
temperature). 

We will say that a spin-flip dynamics is attractive if, for all vertices x of L, the rate 
for the spin flip a x — — 1 — ► a x — +1 is non-decreasing in the number of plus spins in a 
(we will consider only symmetric dynamics, so the roles of plus and minus spins can be 
interchanged). Stochastic Ising models with a ferromagnetic interaction are examples of 
attractive dynamics. A theorem of Harris [221123 states that attractive dynamics preserve 
the FKG property, i.e., if one starts with a measure Pq that satisfies the FKG inequality 
and applies to the spin system an attractive dynamics, the measure Pt describing the spin 
system at time t still satisfies the FKG inequality. 

In particular, this result can be applied to the constrained Glauber dynamics of Sec- 
tion^ to deduce that the limiting (as t — > 00) measure P^ satisfies the FKG inequality. 

2 A Constrained Glauber Dynamics in ID 

We consider the following ferromagnetic Ising chain with single spin-flip (Glauber) dy- 
namics (see [TS] for more details) with (formal) Hamiltonian 



and flipping rates W determined by the energy difference between the configurations after 
and before the proposed move, i.e., 



P(£i\£2)>P(£i). 



(3) 




(4) 



n 



<T n ) = ^SH, 



(5) 



with 



5H = 2(a n _! + a n+l )a n E {-4, 0, 4}. 



(6) 
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Detailed balance with respect to (HJ) at inverse temperature f3 imposes the condition: 



(7) 

Choosing time units so that W_4 = 1, we have W4 = e~ 4/3 . As in (THj, we are interested 
in the zero-temperature case, so that W4 = 0. The rate Wo remains a free parameter, 
which we choose to be zero (corresponding to the zero-temperature limit of the Glauber 
dynamics). Setting Wo = corresponds to allowing only spin flips that lower the energy, 
therefore the only possible moves, happening with rate 1, are flips of plus spins surrounded 
by minus spins or minus spins surrounded by plus spins: 

- + ► , +- + ^ + + +. (8) 

The blocked configurations (i.e., the absorbing states of the dynamics) are those where 
the unsatisfied bonds (i.e., bonds between spins of opposite sign) are isolated. 

We consider the deep- quench situation, where the system is prepared at infinite tem- 
perature and the temperature is then decreased to zero instantaneously. This corresponds 
to an initial configuration chosen randomly from a symmetric Bernoulli product measure, 
i.e., with 

^(O) = +1 with probability 1/2 , . 

<r n (0) = -l with probability 1/2 ^ ' 

where a n (t) is the value of the spin a n at time t. We call P t the distribution of the spin 
configuration cr(t) = {cr n (t)} nG z at time t, and denote by P^ the limiting distribution 
obtained as t — > 00. 



2.1 Checking Edwards' Hypothesis 

In [T5] . Poo is compared to the uniform distribution P un if on blocked configurations, 
corresponding to an ensemble where all blocked configurations have the same weight. The 
conclusions reached there, using exact results on the statistics of the blocked configurations 
reached by the system, invalidate Edwards' hypothesis for this particular model. 

Here, we confirm those results by rigorously proving that P^ 7^ P U nif, but the main 
goal of this section is to introduce, via a simple specific example, a general criterion 
for comparing the limiting distribution P^ of a spin system subjected to an attractive 
dynamics to the uniform distribution P un if on the absorbing states. The general strategy 
is described in Section El where we also give further applications. 

Consider all blocked configurations of the spin chain such that <r± 2 = c±3 = +1- It is 
easy to see that such blocked configurations are of only four different types: 

A a _ 1 = a = a 1 = +l ... + + + + + + + ... 

B cr_i = +1, ct = cri = -1 ... + + + -- + + ... 

C (T_i = (T = -1, 01 = +1 ... + + -- + + + ... 
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D (J_i = (Tq = 0\ = — 1 



+ + + + 



Under the uniform distribution on blocked configuration, the occurrence of each type has 
equal probability; therefore, conditioned on having a± 2 = <r± 3 = +1, 

Pumf^o = +1 | <7 ±2 = a± 3 = +1) = P uni f(A I a± 2 = a± 3 = +1) = 1/4. (10) 

On the other hand, Harris' Theorem (see Section [TT]) applied to this specific dynamics 
implies that P^ satisfies the FKG inequality, so that we have 

P oo (a = +1 | a± 2 = cr ±3 = +1) > P^fo = +1) = 1/2, (11) 

where the equality follows from the ± symmetry of the dynamics and the initial distribu- 
tion. The last two equations show that P^ cannot be the uniform distribution P un if. 

3 The General Strategy 

The strategy we used for the constrained Glauber dynamics of the previous section can be 
generalized to any attractive, symmetric Ising dynamics with locally stable configurations 
(later, we will also give an application where there are no locally stable configurations - see 
Example 4 in Section ETTJ) . with initial configuration chosen from a symmetric distribution 
that satisfies the FKG inequality. For simplicity, we restrict our attention to nearest 
neighbor models; in this context, by the existence of locally stable configurations we mean 
that there are finite subsets G of the lattice L such that, if a x (to) = +1 (—1) Vx G G, 
then a x (t) = +1 (—1) Vx G G, Wt > to. When this is the case, we say that the spins in 
G are stable and we call G a stable set. If G is a smallest set with this property (there 
could be more than one, with different shapes), we call it a minimal stable set. 

Some more notation is needed before we can proceed with the general strategy and 
further applications. Given a subset A of L, we call exterior boundary d e A of A the set of 
vertices x ^ A that are adjacent to a vertex in A, and interior boundary <9jA of A the set 
of vertices x G A that are adjacent to a vertex not in A. 

We are now ready to explain the general strategy; in the next section we will illustrate 
it with some examples. Let G\ and G 2 be two distinct minimal stable sets both containing 
the origin (0 G G\ PI G 2 ) and denote by G = Gi U G 2 their union. Let L be a (finite) 
stable set such that G fl L = and d e G C L (in words, G is "surrounded" by L). G±, G 2 
and L should be chosen so that {G \ G\} U L and {G \ G 2 } U L are stable sets. Notice 
that, since G\ and G 2 are minimal stable sets, G \G± and G \G 2 are smaller than any 
minimal stable set and therefore are not stable sets. 

Now it is easy to convince oneself that, conditioned on the spins in L all being plus, 
there are only four possible types of blocked configurations: 

1. All the spins in G are plus. 

2. All the spins in G are minus. 
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3. The spins in G\ are minus and those in G \ G\ are plus. 

4. The spins in G2 are minus and those in G \ G2 are plus. 

This implies that, conditioned on all the spins in L being plus, the uniform distribution 
on stable configurations assigns probability 1/4 to the event that the spin at the origin is 
plus (corresponding to case 1 above). 

On the other hand, if we consider a symmetric, attractive dynamics with initial con- 
figuration chosen from a symmetric Bernoulli product measure, conditioned on the same 
event (all the spins in L being plus), the limiting distribution P^ must assign probability 
at least 1/2 to the event that the spin at the origin is plus, which clearly shows that P^ 
cannot be the uniform distribution. 

3.1 Higher Dimensional Examples 

Here we present some examples in dimension higher than one where we can use the 
method described above to rule out the uniform distribution. All we have to do is choose 
the sets G%, G 2 , and L appropriately. We will consider zero-temperature dynamics such 
that a spin flips at rate 1 if it disagrees with a strict majority of its neighbors and at 
rate otherwise. As in the example of Section El we will always start with a Bernoulli 
symmetric product measure (see ©), corresponding to the deep-quench situation. We 
note that exact results are usually not available for higher-dimensional models, which 
limits the range of applicability of the methods used in ^J] to one-dimensional models. 

Example 1: Zero-temperature dynamics on the ladder Z x {0, 1}. 

The blocked configurations are such that each spin has at least two neighbors of the same 
sign; squares are minimal stable sets. We choose G\ and G2 to be the sets of vertices of 
the two squares containing the origin {0} x {0} (the shaded squares in Figure [l]) and L 
to be the set of vertices {±2, ±3} x {0, 1}. 

Conditioning on the increasing event £ = {cr±2 = <J±3 = cr' ±2 = <r' ±3 = +1} (see 
Figure HJ), it is easy to see that P un if(cr = +1 | £) = 1/4, because the fact that <r = +1 
implies that er±i = o' ±1 = +1, while there are three possible local blocked configurations 
with cr = —1- On the other hand, P U m/(o"o = +1) = 1/2 by symmetry, so that P un if does 
not satisfy the FKG inequality. 

-3 -2 -1 +1 +2 +3 
j +1+1 +1 + 1 

+1 +1 +1 +1 

Figure 1: The increasing event £ used in Example 1; a n are the spins in the lower row 
Z x {0} and a' n those in the upper one Z x {1}. 
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Example 2: Zero-temperature dynamics on the hexagonal lattice. 

The blocked configurations are again such that each spin has at least two neighbors of the 
same sign; hexagons are minimal stable sets. Let A = GiUG , 2UL = GULbe the set of 
vertices of the portion of hexagonal lattice shown in FigureEl where G\ and Gi are the sets 
of vertices of the two shaded hexagons containing the origin and d e G C L = A \ G = <9;A. 

Conditioning on the increasing event 8 = \o y = +1, G L = <9jA} that all the spins 
in L = <9jA are +1, it is easy to see that Punif{ a o = +1 |^) = 1/4, because the fact that 
cr = +1 implies that a y = +1 for all y G G, while there are three possible local blocked 
configurations with <7o = —1. On the other hand, P U m/(c"o = +1) = 1/2 by symmetry, so 
that Punif does not satisfy the FKG inequality. 




+i +i 

Figure 2: The increasing event £ used in Example 2. As indicated, the vertices of the 
interior boundary <9jA of A all have spin +1. 

Example 3: Zero-temperature dynamics on 7L d . 

For simplicity, we consider the two-dimensional case d = 2, but the same reasoning works 
for all d > 2. In two dimensions the blocked configurations are again such that each 
spin has at least two neighbors of the same sign; squares are minimal stable sets. Let 
A = GiU(j2UL = GULbe the set of vertices of the portion of square lattice shown in 
Figure H2 where G\ and G2 are the sets of vertices of the two shaded squares containing 
the origin and d e G C L = A \ G = <9jA. 

Conditioning on the increasing event 8 = {a y = +1,V?/ G L = <9jA} that all the spins 
in L = diA are +1, it is easy to see that P un j/(cro = +1 | £) = 1/4, because the fact that 
do = +1 implies that o y = +1 for all y G G, while there are three possible local blocked 
configurations with a = — 1- On the other hand, P U nif{&o = +1) = 1/2 by symmetry, so 
that P un if does not satisfy the FKG inequality. 

Example 4: Zero-temperature dynamics on the Cayley tree of degree 3. 
This last example is interesting because, contrary to all the previous ones, there are no 
locally stable configurations (the only stable structures are doubly-infinite plus or minus 
paths). Nonetheless, the criterion described in this paper can still be used. 

With reference to Figure 01 conditioning on the increasing event £ = {a Xl = a X2 = 
a X3 = +1 and xi,X2,x^ belong to doubly-infinite +1 paths that do not contain 2/1,1/2,1/3}, 
it is easy to see that, while P un if(&o = +1) = 1/2 by symmetry, P un if(<7 — +1 1 £) — 1/4, 
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Figure 3: The increasing event S used in Example 3. As indicated, the vertices of the 
interior boundary <9jA of A all have spin +1. 



because if e>"o = +1, then a Vl , a V2 and all forced to be +1. Therefore, once again 

Punif does not satisfy the FKG inequality 




X 3 

Figure 4: A portion of the Cayley tree of degree three. 



4 More on the Constrained Glauber Dynamics 

In [T5] . Poo is compared to two distributions: the uniform distribution P un if on all blocked 
configurations, already discussed in Section 12.11 and an ensemble Pe& where all blocked 
configurations with a given energy density are taken with equal weight. Formally, we can 
write 

P Ed (a) oc e-^\ (12) 

where a is a blocked configuration and E(a) its energy density, and [3 C = 1/T C is the 
inverse of the configurational temperature 0- 

The inverse temperature f3 c that appears in (|T2*|) is itself a function of the energy density 
E, and can be computed using an explicit expression for the configurational entropy 5c 
(see ^3] and references therein) 

S C (E) = E\n {-2E) + In (1 - E) - In (1 + E), (13) 
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from which we obtain 

To derive ()13|). consider a finite chain of iV spins with periodic boundary conditions. 
The blocked configurations such that exactly n bonds are unsatisfied (i.e., are between 
spins of opposite sign) have energy density E(y) = —(N — 2n)/N = —(1 — 2u), where 
v = n/N, and their number is 

Af(N,n)= ( N ~ n )- (15) 

Indeed, this is the number of ways of inserting n unsatisfied bonds between the N — n 
satisfied ones, in such a way that the unsatisfied bonds are isolated. It is now easy to see 
that jt In (J\f(N, n)) — > S c (E(u)) as N, n — > oo and n/N — v is kept constant. 

The translation-ergodicity of the model and the translation-invariance of the energy 
density imply that, with P -P r °bability 1, in the thermodynamic limit, the energy density 
of the final blocked configuration is a deterministic constant, which can be computed by 
solving exactly a dynamical equation for the densities of clusters of consecutive unsatisfied 
bonds. This is done in jT^I (see also (THl E01 EI]), where the energy density for the 
deep-quench situation that we are considering in this paper (see (JHJ)) is shown to be 
E = — 1 + e _1 ~ —0.63212. Therefore, using (|14|). the appropriate value of the inverse 
temperature in (jl2J) for the deep-quench situation is (3 C ~ 0.4895. 

Having an exact calculation for the configurational inverse temperature (3 C , we can use 
our general criterion to check whether = Pecl or not. A simple calculation is sufficient 
to rule out this possibility. Let Pecl{.-) = PEd{- \ c±2 = cr±3 = +1) be the distribution Pe<i 
conditioned on having <j±2 = cr±3 = +1. Then, we have the following straightforward 
relations (with reference to the events A, B, C, and D of Section l2~Tj) 

P Ed (B) =e- 2 ^P Ed (A), (16) 
P Ed (B) = P Ed (C) = P Ed (D), (17) 
P Ed (A) + P Ed (B) + P Ed (C) + P Ed (D) = 1, (18) 

from which it follows that 

P UA) = TT jp Wc - (19) 
Identifying P^ with P Ed would imply, using (fTTj). 

1 + 3e" 2/3c < 2 (20) 

and thus 

Pc > ^inS ~ 0.5493, (21) 

which contradicts the value (3 C ~ 0.4895 corresponding to the deep-quench situation that 
we are considering. 
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